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Using magnetic stripes to stabilize superfluidity in electron-hole double monolayer graphene 
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Experiments have confirmed that double monolayer graphene does not generate finite temperature electron- 
hole superfluidity, because of very strong screening of the pairing attraction. The linear dispersing energy 
bands in monolayer graphene block any attempt to reduce the strength of the screening. We propose a hybrid 
device with two sheets of monolayer graphene in a modulated periodic perpendicular magnetic field. The field 
preserves the isotropic Dirac cones of the original monolayers but reduces the slope of the cones, making the 
monolayer Fermi velocity vf smaller. We demonstrate that with current experimental techniques, the reduction 
in uf can weaken the screening sufficiently to allow electron-hole superfluidity at measurable temperatures. 

PACS numbers: 71.35.-y, 73.22.Gk, 74.78.Fk 


I. INTRODUCTION 

The transition temperatures for electron-hole superfluid¬ 
ity in thin parallel conducting sheets of electrons and holes 
are expected to be high because the electron-hole pairing 
is Coulombic and is strong compared with conventional su¬ 
perconductors. This has led to suggestions of superfluidity 
at room temperatures in double electron-hole monolayers of 
graphen&l, but Re f.|3 showed that strong screening of the pair¬ 
ing attraction in this system tends to suppress finite tempera¬ 
ture superfluidity. Here we propose a hybrid double mono- 
layer graphene device designed to boost the pairing attraction 
by reducing the effects of the screening, and we demonstrate 
that this can lead to magnetically induced superfluidity. We 
use electronic band structure engineering, coupling periodic 
real or pseudo magnetic fields to double electron-hole mono- 
layers of graphene separated by a thin insulating barrier. The 
resulting quantum properties of the device stabilize macro¬ 
scopic quantum coherence and allow, in a solid state electronic 
device, the tuning of the strength of the many-body correla¬ 
tions and the related superfluid properties. Previously, such 
tuning has only been possible in ultra-cold fermionic atoms^. 

Both theory^il^ and experiment^ have established that con¬ 
ventional electron-hole double monolayer graphene does not 
generate finite temperature electron-hole superfluidity be¬ 
cause of strong screening of the electron-hole pairing. This 
originates from the linear Dirac cones of the monolayer 
graphene band structure, e±(k) = ±?iuF|k|, with constant 
Fermi velocity vf, that makes the Fermi energy Ep domi¬ 
nate the average Coulomb interaction (Vcoui)- The resulting 
small interaction strength parameter, = {Vcoui)/Ep = 
e^/ (hvpK) that is fixed independent of density, leads to strong 
screening that makes pairing too weak for finite-temperature 
superfluidity to occur—. Only when = 2.3 does 

the pairing become sufficiently strong for a large superfluid 
gap A to open discontinuously and suppress the screening. 
With graphene on a hexagonal boron-nitride (h-BN) substrate 
of dielectric constant k ~ 3—4 , < 1. Theoretically, weak- 

coupled superfluidity could still occur, but at impractically 
low temperatures. Since it would be destroyed by residual 
disorder— we do not consider it further. Other systems have 


been proposed for observing the superfluidity with an pa¬ 
rameter that can be varied with the density. These include two 
sheets of multilayer graphene which have non-linear dispers¬ 
ing energy bandsi^iS, double quantum wells in GaA a^°’*\ and 
hybrid GaAs-graphene structures'^. 

In this paper we propose use of a periodic magnetic 
field applied perpendicular to double electron-hole monolayer 
graphene in order to reduce the slope of the monolayer Dirac 
cones while preserving their isotrop}^—. This reduces vp 
and increases the value of r^. If Vg can be increased to 
Tg > then finite-temperature superfluidity can occur—. 

The Fermi velocity in monolayer graphene can also be renor¬ 
malized, but non-isotropic ally, by applying a one-dimensional 
potential superlattice in the layer—. 


II. METHODS 
A. Reducing the Fermi velocity 

Consider a magnetic field in the ^-direction, perpendicular 
to the monolayers, as a periodic array in the x-direction of 
rectangular magnetic barriers and wells of height = ±B 
and width dp- The field could be generated with a periodic 
array of ferromagnetic stripes placed on top of the graphene 
monolayers (Fig.[T]i. 

Since the average of the magnetic flux is zero across the 
unit cell of the periodic field, the main effect of the field is 
to modify the monolayer band structure. With zero flux, the 
results are insensitive to fine details of the magnetic profilal^. 
At low energies the de Broglie wavelengths of the quasipar¬ 
ticles are much longer than the length scale of the magnetic 
field variation, so the magnetic profile can be approximated by 
a periodic square wave magnetic fieldl^. We assume smear¬ 
ing of the magnetic barriers is much greater than the lattice 
spacing. The smoothness of the vector potential on a micro¬ 
scopic scale means we can neglect intervalley scattering and 
use single-valley continuum Dirac-Weyl theory. The Zeeman 
effect is very small in graphene, and electron-hole pairing is 
insensitive to relative spin orientation, so we can neglect spin 
effects induced by the magnetic field. 
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FIG. 1: (Color online) Possible realization of the device. An array 
of ferromagnetic stripes with periodicity 2ds placed on top of two 
monolayer sheets of graphene, produces a periodic magnetic field 
B, ~ ±B. 


With this magnetic profile, the vector potential in the Dirac 
Hamiltonian can be fixed by the Landau gauge and chosen pe¬ 
riodic in the x direction. The spectrum can then be obtained 
by a standard transfer matrix approacbi^: Define T (depend¬ 
ing on B, momentum ky and energy level e) as the transfer 
matrix that relates the two-component wavefunction on x to 
its value on x ±2dB- Showing that det T = 1 and introduc¬ 
ing a new momentum due to the periodicity of the superlat¬ 
tice, the eigenvalues of T can be written as so the 

condition determining the band structure is 

Tr[T{ky,e)] = 2cos{2k^dB)- (1) 

Expanding Tr e)] on ky and e, noting that 

Tr[T(A:y,0)] = 2cosh(2fcy(i_B), and expanding cos(2fca;(i_B) 
on kx, we then solve Eq. ([T]i for e. We find that the energy 
dispersion remains linear and isotropic for small momentum, 
but with a reduced velocity aavp < , 

e±(k) = ±fi(adWF)|k|(H-(5(k)) . (2) 

In Eq. (|2l) is a function of c? = dsj^B, the dimensionless 
stripe width, where the magnetic length (.b = yJfic/eB ~ 
26/ \JB[T) nm. An expansion of Eq. ([T]) limits the correction 
term in Eq. I<5(k)| < d^k^/6. The decrease in Eermi 
velocity in the small and large d limits 

ad ~ 1 - d^/60, d < 1 (3) 

2c? j2 ft 

ad ^ —;=e , d » 1 . (4) 

At large densities, Ep eventually passes out of the first 
energy band of the periodic magnetic field into the band 
gap where the linear spectrum approximation is no longer 
valid. The energy width Eb of the band can be numer¬ 
ically calculated by solving Eq. O at the boundaries of 
the Brillouin zone, kx = ±7r/2dB, i.e. Tr[T(0, = 

-2, with Eo = hyp/lB ~ 22^5 (Tesla) meV. We 
find Egad/d < Eb < i5o7rad/2d. The lower limit is 
valid for large d. It corresponds to the spectrum along 
kx as e{kx) = h{adVF)y/[l - cos(2A:a:dB)]/2/dB ^ 


h{adVF)kx{i — d^gk'^/Q), the largest deviation from linear¬ 
ity. The upper limit is valid for small d. It corresponds to 
a completely linear spectrum, e{kx) = hadVFkx- Remark¬ 
ably, along the ky direction the spectrum is found to be al¬ 
most linear for all e < Eb and for all parameters we use. 
Within linear approximation, we must restrict our results to 
values of Ef < Eb- Even in the least favourable case, when 
the extreme value Ep = Eb is reached at kx = l/ds (in 
the linear approximation), the true spectrum remains close to 
the linear spectrum, e(^) = y/(l — cos 2) /2Ef — 0.84£'f. 
This gives an estimate of the maximum error along kx- In 
the superfluid state, i5(k) must also be small compared with 
the energy gap Amax since the gap excludes single-particle 
states lying less than A^ax above Ep, or equivalently, using 
the limiting value at the Brillouin zone edge, |i5(k)| < 7r^/24, 

7r^ad/48d < A^ax/^^o • (5) 

Eor d 1, Eq. 0 is always satisfied since a^ ~ de 


B. Enhancement of superfluidity by magnetic field 


The effective Hamiltonian for the two monolayer sheets of 
graphene in the presence of the periodic perpendicular mag¬ 
netic field is, 


k7 


qk7 


7t 7t 

7+a^-k+a'-k'-i-a'"-k'+a 


H -(6) 


The single-particle energy bands for the modified Dirac cones 
of the conduction band (electrons) and valence band (holes), 
^k “ TCtd't’F |k| — fJ-, are measured from their respective chem¬ 
ical potentials ±/i, where 7 = 1(—1) labels the electron (hole) 
sheet. The and are creation and destruction operators 
for electrons and holes. Spin indices are implicit. is the 
screened electron-hole interaction. 

The mean-field equations at zero temperature for the 
momentum-dependent superfluid gap functions A^, and for 
equal electron and hole densities n+ = n- = n are, 



E'^ = Y^k ^ + ^k 9s{gv) = 2 are the spin (pseudospin) 
factors, and H is the sheet area. We retain only the s-wave har¬ 
monic in the graphene form factor F/*/) =1/2 which comes 
from the overlap of the single-particle wave-functions in the 
strong-coupled regime^. 

We self-consistently calculate the screened electron-hole 
interaction within the random phase approximation 

(RPA) in the zero temperature superfluid state ^A^° . The most 
favourable conditions for pairing are at small interlayer sepa¬ 
rations D on the scales of both the effective Bohr radius and 
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FIG. 2: (Color online) Interaction parameter Va for magnetic field of 
periodicity 2d for dielectric constant ac = 3 (monolayers embedded 
in a h-BN substrate and k = 2 (free standing monolayers, separated 
by h-BN). Superfluidity occurs for > rr“ = 2.35^. 


FIG. 3: (Color online) Maximum superfluid gap Amax for different 
values of the magnetic field B, as a function of sheet densities n. 
Sheet separation is D = 2 nm. At arrow points: Amax in Kelvin. 


the inverse Fermi momentum in each layer. In this case 
qD 1 and the interaction reduces to. 


yeh 

q 


Vne 


-qD 


1 + 2vqIV{q) 


(9) 


Vq = —2'Ke^l{Kq) is the unscreened Coulomb interaction. 
U{q) = ~ is the 

sum of the normal (intralayer) and anomalous (interlayer) po¬ 
larizabilities for the superfluid state calculated with the lin¬ 
ear energy spectrum. The gap equation [Eq. (|7]i] is inde¬ 
pendent of density when expressed in units of Ep and kp, 
with the exception of the e~‘^^ factor in Eq. 0 for 
With increasing density, this factor weakens Rescaling 
q = adq gives n(g, aavp) = aj^n(g, vp), and, for qD <C 1, 
V^^{adVp, n) = a'^Vq’^{vp,adK). Equations Q and |8] thus 
remain the same, but with (adn) replacing k. The interaction 
parameterr^, = [l/ad\[e^ / {finvp)] increases by a factor 


III. RESULTS 


Eigure |2] shows that by tuning ad with d using a periodic 
magnetic held, Vg can be increased above the value r™*®* 
needed for superfluidity, k = 3 corresponds to monolayers 
embedded in a h-BN substrate, k = 2 corresponds to a free 
standing system with the two monolayers separated by h-BN. 
Eor d < dmin ~ 2Y^logK + l, Tg < r™*®', and the superfluidity 
is killed by strong screening. Eor k = 3(2), d^m ~ 3.1(2.7). 

The renormalization of the band structure is the main effect 
that drives the superfluidity in the double monolayer system. 
We And whenever superfluidity occurs, the electron-hole pairs 
of the superfluid ground state are compact compared with their 
spacing, making them approximately neutral, and the elec¬ 
trons and holes have opposite wave-vectors. Thus screening 
effects and effects of the magnetic held on the orbital degrees 
of freedom should be small compared with the primary effect, 
the renormalization of the Eermi velocity. 


Eigure 0 shows the maximum superfluid energy gap Amax 
at zero temperature for different values of the magnetic held 
B, as a function of sheet density n, for two monolayers sepa¬ 
rated by D = 2 nm and embedded in a h-BN substrate. The 
gaps Amax are of the order of several hundred Kelvin. Amax 
decreases with increasing n, due to the e~‘^^ factor in 
[Eq.ig. This eventually results in no solution to the gap equa¬ 
tion. However, we terminate the curves in Eig. 0when Ep 
reaches Ep, and this occurs before such a density is reached. 
Ep = Ep thus gives a lower limit on the maximum density 
for the superfluid phase, n ~ 5 x B{E)/d? cm“^. 

Eor all the gaps shown in Eig.0 Amax ^ Ep, leading to 
a strong suppression of the screening. At higher densities, the 
very strong screening would kill the superfluidity before the 
system can enter the BCS regime^. However, with increas¬ 
ing density, and well before Amax can drop to Ep, the Eermi 
energy reaches the band edge Ep, where the curve must be 
truncated. 

To maximize the density range for superfluidity, the mag¬ 
netic held B should be made large. In the ferromagnetic 
stripes shown in Eig. [T] the maximum magnetic held is i? ~ 
1.2 T, with magnetic length ip 2A nm. Eor k = 3, this 
corresponds to a minimum stripe width dp = dminis = 74 
nm, readily attainable experimentally. However for i? ^1.2 
T, the magnetic band width is narrow so that Ep reaches Ep 
at low densities (see Eig.0. 

To obtain superfluidity at higher densities, a deforma¬ 
tion of the graphene layer can be used to produce much 
larger pseudo-magnetic fields. Periodic deformations of 
the graphene layers generate a fictitious vector potential 
which can produce a periodic pseudo-magnetic held in the 
layers^-^. 

The pseudo magnetic held induced by the strain changes 
sign for the two valleys. The Dirac cones at the Brillouin 
zone points K and K' will experience an alternating mag¬ 
netic held in both cases, but with a tt shift in phase. However, 
in the absence of intervalley scattering, as considered in this 
work, a global tt shift in the zero-flux magnetic field does not 
affect the renormalization of the Eermi velocity of the two 
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through surface acoustic waves (SAWs)21, using two inter¬ 
digital transducers and a piezo-electric substrate. For wave¬ 
lengths A ~ 80 nm and typical piezo-electric materials used 
in SAW devices, the frequency ~ 50 GHz. 

The Fourier transform of the interaction in Eq. (|9]l in the 
presence of graphene corrugations should be evaluated in the 
periodic curved geometry. However, in the strong coupling 
regime, the electron-hole pairs have dimensions along the 
layers comparable to the deformation wave-lengths, so the 
curved geometry corrections to the interaction will be small. 

A. Transition temperature 


FIG. 4: (Color online) (a) Periodic deformation profile /i of a 
graphene sheet, (b) Induced strain in the sheet, (c) Induced pseudo- 
magnetic field in the sheet. 


Dirac cones. Therefore, the evaluation procedure for deter¬ 
mining the RPA screening and the density of carriers, using 
two equivalent renormalized Dirac cones, is the same as for 
the ferromagnetic stripes. 

Let us consider a periodic modulation of the graphene layer 
with height profile [Fig. 15 a)] along the zig-zag direction of 
the lattice, h{x,y) = ho cos (27ra;/A). ho and A are the 
modulation amplitude and wavelength. The pseudo-potential 
A = {Ax, Ay) at lattice point a is defined as; 

AZ+%A’^ = — V . (10) 

" evp , 

I3^nn{a) 


The K point of the Brillouin zone is at K, Rap is the distance 
between atoms a and /3, and the hopping amplitude tap cou¬ 
ples the Pz orbitals on neighboring atoms. The corresponding 
pseudo-magnetic field is 5 = V x A. 

Figure |5a) shows the profile of the periodic deformation 
along the zig-zag direction for A = 80 nm and amplitudes 
ho. Figure |5b) shows the induced strain, and Fig. |5c) the 
induced periodic pseudo-magnetic field. The periodicity of 
the strain profile and the pseudo-magnetic fields is one-half of 
A, so A = 80 nm corresponds to a stripe width hs = 20 nm. 

To explore the modification of the electronic properties in¬ 
duced by an out-of-plane deformation, we have calculated the 
local density of states (LDOS) using the tight-binding Hamil¬ 
tonian with a spatially varying hopping amplitude induced by 
a spatially varying inter-carbon distance. The LDOS is cal¬ 
culated through a Chebyshev expansion of the single particle 
Green’s functio n^^i^^ . Lor A = 80 nm, we find that for de¬ 
formation amplitudes up to ho < 4 nm, we recover the same 
phenomenon in the energy dispersion as that for the real mag¬ 
netic field, an isotropic linear dispersion with the slower Lermi 
velocity adVp- Ligure|5c) shows that an amplitude hp ~ 2.6 
nm generates the large pseudo-magnetic field i? ~ 20 T. This 
leads to a much wider magnetic band width Eb than is possi¬ 
ble for ferromagnetic stripes, and Ep does not reach Eb until 
densities n > 10^^ cm“^. An experimental realization would 
be to deposit graphene on a substrate that can be strained 


Over the range of parameters we are considering, we find 
that the electron-hole superfluidity is always in the strong¬ 
coupling regime, well inside the crossover regime of the BCS- 
BEC crossover. The transition temperature calculated within 
the mean held approach, Tm /, will be much larger than the ac¬ 
tual transition temperature for the onset of phase coherence. A 
lower bound on the transition temperature in two-dimensions 
is given by the Kosterlitz-Thouless (KT) temperature, 

Tkt = {t^/2)ps{Tkt) , (11) 

where Ps{T) is the superfluid stiffness. Ps{T = 0) = Ep /Ap 
and, when kpD is small, Ps{T) falls off slowly with T for 
T Ainax- Hence taking Ps{T) ~ Ps(0), we obtain Tkt = 
Ep/9> for Ep/d> Amax- 

Since in general Tkt ^ Tmf, the mean field superfluid 
gap will be insensitive to T for T < Tkt- By a similar argu¬ 
ment, the RPA screening polarization bubbles are only weakly 
affected by finite T < Tkt- Thus, in the superfluid calcula¬ 
tions we can take T = 0. 

Lor superfluidity to occur, we recall that the array spacing 
d ^ 3 (Lig. 15. Lor these values of d, the linearized Eq. (|5 
is valid whenever Ep < Eb is satisfied. This inequality 
establishes an upper bound on the maximum transition tem¬ 
perature, Tkt < T'jJ” — 35 (Tesla) exp {—d'^/A) K, the 
equality occurring for Ep = Eb- Ligure|5] shows this maxi¬ 
mum as a function of both B and d. 




FIG. 5: (Color online) (a) as a function of B. An'ay spacings 
d = 3.1 and 2.7 give maximum transition temperature for embedded 
and free-standing system, respectively, (b) as a function of d. 
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IV. CONCLUSIONS 

We have proposed an electron-hole double monolayer 
graphene system, in which a quantum phase transition to a 
superfluid is induced by a periodic magnetic field. Electron- 
hole pairing in conventional double monolayers is known to 
be severely weakened by screening within the layers which 
kills the superfluidity. If the pairing can be made strong, then 
a large superfluid gap is known to open up, destroying the 
low-energy single-particle excitations that cause the screen¬ 
ing. We show how a periodic magnetic field applied perpen¬ 
dicular to the monolayers could be used for this purpose; The 
field preserves the isotropic Dirac cones of the original mono- 
layers but reduces the Fermi velocity in a tunable way, shift¬ 


ing the system parameters into the strongly coupled pairing 
regime where screening is sufficiently weakened for finite- 
temperature electron-hole superfluidity to occur. 
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